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Abstract 

We contribute to the method of trigonometric series for solving 
differential equations of certain non linear oscillators. 
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1 Introduction 

The non linear nonharmonic motion of an oscillator may be given by the 
following differential equation 

u" + u?u = -(3u 2 (1) 

(3 and u being constants, with initial conditions 

u(0) = oo, u'(0) = (2) 

To solve this problem A. Shidfar and A. Sadeghi [1], have given two 
series solutions. They describe a general approach in which the differential 
equation, rather than the solutions series , is majorized. 
Notice that if we write ao = — Hj- then 

u(t) = -j (3) 
is a trivial solution of (1) and (2). 

They gave a series solutions of (1) and (2), which includes (3) as a special 
1 2000 Mathematical Subjects Classification, 34A20. 
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case. 



By writing 

co 2 

U = V ~2P 

the problem becomes 

v" + 0v 2 = ^ (4) 

under the initial conditions 

f u(0) = a + g 

\ w'(0) =0 V ; 

The method of [1] consists to solve equations (4) and (5) in the form 

v(t) = Co + c\ sin ut + C2 sin 2 ut + C3 sin 3 + (6) 

where q, z = 0, 1, 2, .... are coefficients to be determined by the substitution 
of (6) in (4). 

In fact, to = y where T is the period of the solution, which can be expressed 
in terms of the Weierstrass function p(z, 2T, 2T'). 
So, we find that 

2uj 2 c 2 + 0c 2 o = —. 
For n > 1, the recursion formula for these coefficients are 

n 

(n + l)(n + 2)c n+ 2 = n 2 c n ^ C r c n-r- (7) 

U r=0 

Equations (5) and (6) imply that c\ = 0. Relations (7) yields 

C3 = 0, c 5 = 0, 

The even order coefficients simply are 

CO 2 

co = a + ^, 

C2 = -^(oo 2 + a o 0), 

( 2 , flW 3 ao/3s 
C4 = -6^ O0( ^ +ao/?)( 4-2^ } ' 

/ 2 , fl \/3 ao/? v 2ao/? A /3a 2 , , 2 , fl \2 
C6 = -180^ ao( " + ^ )( 4"y )(15 -^ ) -i2P (w +ao/3) ' 
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etc. 

The coefficient cq follows from the condition (5). The solution for the 
equations (4) and (5) can now be written as 

u(t) = a - T^r + a o/3) sin 2 ut + .... (8) 
Relations of the coefficients and further induction show that C2i, % = 1, 2, ... 

2 

all vanish for cto = So, the trivial solution (3) is included in (6) as a 

special case. 

2 Convergence of the solutions 

We now show the convergence of these series. In [1] one proved the following 



Lemma 1 The serie (6) solution of Equation (4)-(5) is absolutely 
convergent for all t. 

Proof We firstly note that if c$ > 0, c 2 > and (3 < 0, then all 
coefficients c n in the serie expansion (6) are positive. Indeed, we may write 

^((n + l)(n + 2)c„ +2 = 7)2 c « ~ ^(^Z c «) 2 + ^ c o> 

n>0 n>0 W n>0 W 

or 

-/3(^ Cri ) 2 + co 2 ^nc ra = -/?cg-2u; 2 C2 . (9) 

n>0 n>0 

Since the right hand side of (9) is finite and Cj are positive, the series 
J2n>o c n converges. Following [1], if we put 

c o =1 c o I) A = 0, c' 2 =| c 2 | 

and for n > 2 

° n+2 = (n+l)(n + 2) Cn + w 2 (n + l)(n + 2) ^ CrCn " r ' 

then the series X)n>o c n converges. Since | c n |> c^, it follows that the 
solution series (6) is absolutely convergent, and hence the series expansion 
solution of (l)-(2) converges for all t. 

We notice that we may deduce Lemma 1 a previous result concerning 
Equation (4). 

We have shown that the coefficients verify a more general properties. Indeed, 
we have [4] 
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Lemma 2 For any positive number e small enough (but e / 0), 
there exists a positive constant k verifying 

k < 4L 4 -i 

LO 2 4 

such that the coefficients c n of the series expansion (6) solution of the dif- 
ferential equation (4)- (5) satisfy the inequality 

I C n \< (10) 
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Proof We first notice that Lemma 2 gives an optimal result, because 
our method do not run for e = 0. 

The coefficients c n of the power series solution, satisfy the recursion formula 
(7). We shall prove there exist two positive constants k > 0, and a > 1, 
such that the following inequality holds 

k 

\ c n < — 

for any integer n > 1. Suppose for any n < p , we get | c n \< ■ n %. In 
particular, it implies that 

k 2 k 2 
V r r < V - < 



r a (r — v) a ~ (v — l)^ 1 

0<r<p 0<r<p v y > V > 



Equality (7) gives 



P 2 ~ 2^cp (3 ^ 

Cp+2 " (p + l)(p + 2) Cp u> 2 (p + l)(p + 2) ^ ° rCp - 



Thus, if we prove the following inequality 
p 2 - 2-4c k (3 



< 7 KT- ( n ) 



(p + l)(p + 2) p Q u 2 (p + l)(p + 2) (p - I)"" 1 - (p + 2) Q 
so 

I ^2 l< (12) 

Notice that (11) implies that 

p 2 -2^c k < k 



(p + l)(p + 2)p« - (p + 2)"' 



4 



Thus, a necessary condition to (11) holds is : a < |. Inequality (11) is 
equivalent to 

k < ^Pf(p)g(p) (13) 

where 

p p + 2 

= 1 ¥TW« 

By using MAPLE, we are able to prove that f(p) is an increasing positive 
function in p . Moreover, for any p > 1, f(p) is minorated 

Hp) > (In 1 --. 

The function g(p) is such that 

(p 2 - Jrcox^r- 1 



(p+1) 



is a strictly decreasing and bounded function . 
More exactly, we may calculate the lower bound 

. . (3 - 2a) 
P 

Thus, if (3 — 2a) = e > 0, it suffices to choice 

fc <(3 )4 i- Q(3 _ 2a) 

to inequality (13) holds. 



Remark for the case e = : 

Notice that the choice of k depends on a value. 

For a = |, we then prove by MAPLE that the function 

P9(P)=P ^ 

is positive and strictly decreasing to 0. While p 2 g(p) is a bounded function. 
Moreover, it appears that pf(p)g(p) is a decreasing function which tends to 
when p tends to infinity. Thus, our method falls since it do not permit to 
determine a non negative constant k. 



5 



Following [1], it is interesting to write the power series solution for the 

system (4)-(5), 

oo 

v(x) = ]T b n x n . (14) 

n=0 

hp+i = p = 0,1,2,..., 



We find again that 
while 



J 1 



bo = a + -, 
2b 2 + (3b 2 = ^, 

r=n 

(n + 2)(n + l)6 n+2 = -/3 ^ Mn-r, 

r=0 

where n is even and non zero. 

2 

The coefficients b2 P ,p = 0, 1, 2, again vanish for ao = — tj-. We may verify 
that the solutions (6) and (14) are identical. 

The latter method permits to compare approximate solutions of the anhar- 
monic motion of the oscillator. 



3 Another differential equation 

We now examine the following differential equation 

u" + u?u = -f3u 3 (15) 
j3 and u being constants, with initial conditions 

u(0) = ao, «'(0) = 0. (16) 

We put v = and t = ujx. We then obtain from (15) and (16) 

d 2 v 

_ + t , + /3 w 3 = 0, u(0) = l, t/(0) = (17) 

where 3 = 

^ to 

A. Shidfar and A. Sadeghi [2] solved (17) by series method in Sinus power 

v(t) = Co + c\ ■smujt + C2 sin 2 ut + C3 sin 3 ut + (18) 

Here, u = ^ where T is the period of the solution, which can be expressed 

in terms of the Jacobi function sn(z, 2T, 2T"). 

So, 

c = a . 
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For n > 1, we get the recursion formula 

P n n—r 

(n + l)(n + 2)c n+2 = n 2 c„ 2 X X c m c r c„_ m _ r . (19) 

r=0 m=0 

Under some conditions, they proved estimates of the coefficients 

I c n |< R n 
where j=r is a radius of convergence. 

In fact, we may prove an analog of Lemma 2 for this equation. Indeed, we 
have 

n n—r n n—2 n—r 

^ ' ^ ' c m c r c n — m — r = 2cgC n +2coCic n _i+co ^ ' c n c n — m -\- 'y ' ^ ' c m c r c n _ m _ r . 

r=0 m=0 m=0 r=2 m=0 

Then, 

n—2 n—r n—2 n— r— 1 

^ ' ^ ' c m c r c n _ m _ r = ^ ] c r [2coc n _ r + ^ ' c m c n _ m _ r ] 

r=2 m=0 r=2 m=l 



n-2 fc2 
< I °r I I 2 I c 0Cn-r | +7 ~— J 

t^L (n-r- l) a 1 



r=2 

n-2 -, n-2 



< 2c k 2 V — ^ \-k 3 V ; - r < . r + - ^- r 



r=2 v > r=2 

Finally, 



Cp+2 |< 



{p + 2) a 

as soon as the non negative constant k satisfies the inequality 

k + 2c k 2 + 2c ci/c + 3c k 2 + fc 3 < (n + l)fc 



n a " 2 n a (n-l) a (n-2)"- 1 (n-2)"- 1 (n + 2)"- 1 ' 
So, 

1 fc 2 + 3c k 2c + 2c ci (n + 1) 
n°- 2 + (n - 2)"- 1 + (n-l) Q < (n + 2)"- 1 ' 
By using MAPLE, we verify it is possible to find a such constant. Moreover, 
we find again the necessary condition : a < §, since we get w J_^ < (J+2)"- 1 • 

General remarks : It is wellknown from the theory of elliptic 
functions that solutions of equations (4) and (17) are related. This allows to 
express the series expansion of the solution of (17) from a series expansion 
of a solution of (4) and conversely. Indeed, one has 

S 2 

p(z) =C 27FT' 
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where p(z) is the elliptic Weierstrass function and sn(u) is the Jacobi fonc- 
tion , 5 is a constant , only dependent on the initial parameters. 
Notice that series expansion of the sn(u) in sinus power was given in a 
previous paper (see Proposition (2.1) in [3]). 
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